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l/^ , ABSTRACT 

■ The bright southern binary star f3 Centauri (HR5267) has been observed with the 
^-H ' Sydney University SteUar Interferometer (SUSI) and spectroscopically with the ESO 

, CAT and Swiss Euler telescopes at La Silla. The interferometric observations have 

■ confirmed the binary nature of the primary component and have enabled the de- 
termination of the orbital parameters of the system. At the observing wavelength 

. of 442 nm the two components of the primary system have a magnitude difference 

pi i' of 0.15±0.02. The combination of interferometric and spectroscopic data gives the 

following results: orbital period 357.00 ± 0.07 days, semi-major axis 25.30 ± 0.19 mas, 
inclination 67. 4±0. 3 degrees, eccentricity 0.821±0.003, distance 102.3±1.7pc, primary 
and secondary masses Mi — M2 — 9.1 ± 0.3 Mq and absolute visual magnitudes of 
^ i the primary and secondary Miv — —3.85 ± 0.05 and M2V = —3.70 ± 0.05. The high 

' accuracy of the results offers a fruitful starting point for future asteroseismic modelling 

. 5^ ' of the pulsating binary components. 

, Key words: binaries: general — stars: individual (/? Cen, HR 5267) — stars: dis- 

' tances — stars: fundamental parameters — techniques: interferometric — techniques: 

spectroscopic 
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1 INTRODUCTION 

The bright southern star /3 Centauri (HR 5267) is classed 
in the Bright Star Catalogue JHotHeit 198211 as spectral type 
Bl III, and as a visual double of separation I'.'S with a mag- 
nitude difference of 3.2. In 1963, during the early commis- 
sioning observations with the Narrabri Stellar Intensity In- 
terferometer (NSII) jllanbury Brown, Davis fc Allen 1967| l, 
it was found that the instrument's response to /3 Cen was ap- 
proximately half that expected from a single star. From this 
observation it was deduced that the primary of the visual 
double was in fact a binary system with two roughly equal 

* E-mail: davis@physics.usyd.edu.au 

f Present address: Astrophysics Group, Cavendish Laboratory, 
Madingley Road, Cambridge CBS OHE, UK 

I Present address: Jet Propulsion Laboratory, California Institute 
of Technology, Pasadena, CA 91109, USA 

§ CSIRO Division of Telecommunications Sz Industrial Physics, 
West Lindfield, NSW 2070, Australia 



components ( |Hanbury Brown, Davis fc Allen 1974| l. It was 
not possible to determine additional information about (3 
Cen with the NSII, but support for the binary nature of the 
primary came from the spectroscopic observations of Breger 
(1967 ), Shobbrook and Robertson (1968), and Lomb ^^7^. 
They found that the primary probably contained a /3 Cephei 
variable with period of either 0.135 d or 0.157d, and that 
residual velocities for the lines are consistent with a binary 
which Shobbrook and Robertson suggested had a period of 
352 d. The spectroscopic observations led to /3 Cen being 
classified as a single-lined spectroscopic binary but Robert- 
son et al. |1999> have used aperture masking interferometry 
and high resolution spectroscopy to confirm that the pri- 
mary is indeed a binary, to resolve the angular separation 
of the two components, and to establish it as a double-lined 
spectroscopic binary. Ausseloos et al. 12U02II have made a 
detailed spectroscopic study of /3 Cen and shown that it is 
an eccentric binary with two /3 Cep-type components and 
determined the orbital parameters of the system. 



2 J. Davis et al. 



In this paper we report observations of the /3 Cen sys- 
tem made with the Sydney University Stellar Interferometer 
(SUSI) ( Dav is et al. 1999a)l which not only confirm the de- 
tection of the two components of the primary reported by 
Robertson et al. but also enable the orbital parameters of 
the system to be determined. In this paper '/3 Cen' will be 
taken to identify the spectroscopic binary alone unless the 
fainter visual companion is specifically mentioned. We also 
report an update of the spectroscopic study of (3 Cen by 
Ausseloos et al. II2002II and combine the interferometric and 
spectroscopic results to determine the distance to the system 
and the masses and absolute magnitudes of the components. 



2 INTERFEROMETRIC OBSERVATIONS 

SUSI is a long baseline optical stellar interferometer with 
a North-South baseline array. It measures the square of 
the fringe visibility which we term the 'correlation' C 
ijDavis et al. 1999b|l . The correlation for a binary star is 
given by 



C = 



[r? + /32ri + 2/?|ri||r2|cos(a)] 



(1) 



where /3 is the brightness ratio of the two components (de- 
fined so that /9 < 1) and 



2Ji{neuDib/Xo) 

n6uDib/Xo 
2Ji{-KeuD2b/Xo) 

2nbp cos tp 



(2) 
(3) 
(4) 



The symbols in equations (|5J and ^ are defined as 
follows; b is the length of the interferometer baseline pro- 
jected onto the sky, Ao is the observing wavelength, 9udi 
and OuD2 are the equivalent uniform-disk angular diameters 
of the primary and secondary components, and Ji{x) is a 
Bessel function. In equation (|1J p is the angular separation 
of the two components and tp is the angle between the line 
joining the two stars and the projection of the baseline onto 
the sky. 

The rotation of the Earth causes both b and tp to vary 
with time. ip will also vary due to orbital motion of the bi- 
nary system. The result of these variations is to cause the 
cosine term in equation Q, and hence the observed correla- 
tion, to vary with hour angle. The timescale of the variation 
in correlation depends on the baseline employed and, for 
each observing night, the baseline was chosen to give at least 
one, and preferably two minima in the observed correlation. 

In order to obtain accurate measurements of the corre- 
lation it it essential to ensure that any difference in optical 
path length AOPL between the two arms of the interferome- 
ter is small compared with the coherence length of the inter- 
ferometer. In SUSI this is done by means of an optical path 
length compensator (OPLC) IIDavis et al. 1 999al . Ideally an 
automatic fringe tracker would be used to control the OPLC 
and maintain path equality but, at the time of the observa- 
tions reported here, this was not available and the following 
technique was employed. The correlation was measured for 
several values of the optical path difi'erence between the two 



arms of the interferometer bracketing the matched path po- 
sition. The correlation at AOPL = was found by fitting 
a 'delay' curve to the data IIDavis et al. 1999bll . A series of 
delay curves was obtained for each night to give a time series 
of correlation values for AOPL — 0. 

The usual mode of operation for interferometric obser- 
vations is to interleave observations of the programme star 
or binary system with observations of calibrators, stars that 
are either unresolved or of known angular diameter. The 
limiting magnitude of SUSI at the time of the observations 
was B ~ 2.5 and there were no suitable calibrators in the 
vicinity of /3 Cen. The observed correlation from /3 Cen was 
generally changing rapidly due to the cosine term in equa- 
tion (0 and it was desirable to monitor it as continuously as 
possible. In the absence of suitable calibrators and the lack 
of a fringe-tracking system it was decided to omit observa- 
tions of calibration sources and, instead of endeavouring to 
measure the angular sizes of the components as well, to con- 
centrate on the determination of the orbital elements of the 
system. This omission adds some complication to the anal- 
ysis but has minimal implications for the determination of 
the separations and position angles of the component stars. 
With these considerations in mind, the observations were 
made at a wavelength of 442 nm and with baseline lengths 
of 5 m, 10 m, 15 m, 20 m or 40 m. At these baselines the two 
component stars of the system are only partially resolved 
and the correlation for either of the components alone is es- 
timated to be > 0.95 for all the baselines except 40 m where 
the correlation is estimated to be >0.8 (see Section I^J. 

A total of 43 nights of observations of /3 Cen were made 
with SUSI during the period April 1997 to April 2002. The 
dates of the observations and the baselines used are listed 
in Table El 



3 ANALYSIS OF THE INTERFEROMETRIC 
OBSERVATIONS 

The measured values of correlation are not the true values of 
correlation due to instrumental and seeing losses. Temporal 
seeing effects are corrected on-line by the SUSI data acqui- 
sition system jPavis fc Tango 19'96llDavis et al. 1999bt but 
spatial seeing effects and instrumental losses could not be 
corrected in the normal way in the absence of calibration 
source observations. However, since the aim is to determine 
the separation and position angle of the binary, the absolute 
values of the correlation are not important. Even for the de- 
termination of the brightness ratio /3, it is only the ratio of 
maximum to minimum correlation that must be determined 
and this is discussed in Section f3. 21 

A least squares fit of equation Q with Fi = r2 = 7 
was made to the data for each night with 7, p and 9 as 
free parameters where 9 is the position angle of the binary 
system (measured east from north) . The justification for the 
assumption of Fi = r2 or, in other words 9udi = 9ud2, was 
initially based on the finding of Hanbury Brown, Davis & 
Allen 119741 that the two component stars of /3 Cen are of 
comparable brightness. Following the determination of the 
brightness ratio (3 (Section 13.211 the angular diameters of 
the primary and secondary were estimated to be 0.77 mas 
and 0.72 mas respectively based on the angular diameters of 
early- type stars by Hanbury Brown, Davis & Allen 1197411 . 
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Table 1. Measured angular separations (p) and position angles (d) for /3 Cen. Ap 
and A6 are the differences between the measured values and the corresponding p and 
d values for the fitted orbit. The weights are those assigned for the orbital fitting 
program (see Section |5J. Further details are given in the text. 



Date Mean B p Ap 6» A6» Weight 



(UT) 


MJD 


(m) 




(mas 






(mas) 




(deg) 






(deg) 




27/4/97 


50565.66 


10 


21, 


.69 


± 


0, 


,20 


-0.05 


118, 


,07 


± 


0, 


.16 


0.14 


3 


8/5/97 


50576.56 


5 


24, 


.83 


± 


0, 


,39 


0.48 


120, 


,54 


± 


0, 


,52 


-0.24 


3 


9/5/97 


50577.59 


10 


24, 


.03 


± 


0, 


19 


-0.53 


121, 


.52 


± 


0, 


.34 


0.49 


4 


10/5/97 


50578.55 


5 


25, 


.94 


± 


0, 


,39 


1.19 


119, 


.73 


± 


0, 


.42 


-1.52 


4 


29/5/98 


50962.57 


20 


28, 


.07 


± 


0, 


,05 


-0.37 


126, 


.65 


± 


0, 


.05 


0.14 


5 


30/5/98 


50963.53 


10 


28, 


.00 


± 


0, 


,06 


-0.52 


127, 


,45 


± 


0, 


,11 


0.77 


2 


8/6/98 


50972.50 


5 


28, 


.94 


± 


0, 


,35 


-0.24 


128, 


,11 




0, 


,23 


-0.07 


4 


9/6/98 


50973.50 


10 


28, 


.82 


± 


0, 


,06 


-0.42 


128, 


.53 


± 


0, 


.07 


0.19 


5 


23/7/98 


51017.45 


10 


29, 


.68 


± 


0, 


,25 


-0.14 


135, 


.33 




0, 


.27 


0.09 


4 


20/3/99 


51257.64 


20 


12, 


.63 




0. 


,06 


-0.15 


107, 


,54 




0, 


,15 


-0.04 


5 


25/3/99 


51262.63 


20 


15, 


.79 




0. 


11 


0.30 


109, 


,73 


± 


0, 


,25 


-1.35 


4 


26/3/99 


51263.62 


20 


15, 


.81 




0. 


,05 


-0.15 


111, 


,52 




0, 


,10 


-0.12 


5 


27/3/99 


51264.52 


20 


16, 


.37 


± 


0, 


,25 


0.00 


111, 


.59 


± 


0, 


.61 


-0.52 


3 


29/3/99 


51266.61 


20 


17, 


.16 


± 


0. 


,06 


-0.12 


112, 


.88 


± 


0, 


.10 


-0.26 


5 


6/4/99 


51274.67 


20 


20, 


.06 


± 


0, 


,04 


-0.19 


116, 


,27 


± 


0, 


.08 


-0.07 


5 


10/4/99 


51278.75 


15 


21, 


.08 


± 


0, 


.10 


-0.41 


117, 


.45 


± 


0, 


.05 


-0.20 


4 


11/4/99 


51279.62 


15 


21, 


.57 




0, 


,04 


-0.17 


117, 


.84 


± 


0, 


.08 


-0.07 


5 


12/4/99 


51280.60 


15 


21, 


.75 




0, 


,03 


-0.25 


118, 


.14 


± 


0, 


.05 


-0.05 


5 


18/4/99 


51286.64 


15 


23, 


.30 




0, 


,05 


-0.21 


119, 


,84 




0, 


.05 


0.02 


5 


19/4/99 


51287.58 


10 


23, 


.18 




0, 


,12 


-0.53 


120, 


,39 




0, 


.14 


0.34 


4 


22/4/99 


51290.59 


10 


24, 


.38 




0, 


,13 


0.02 


120, 


.66 




0, 


.14 


-0.12 


4 


23/4/99 


51291.59 


5 


25, 


.09 


± 


0, 


.46 


0.53 


120, 


.47 




0, 


.50 


-0.54 


4 


28/4/99 


51296.60 


10 


25, 


.52 


± 


0, 


,20 


0.03 


121, 


,98 




0, 


.24 


-0.15 


4 


29/4/99 


51297.60 


5 


26, 


.57 


± 


0, 


.37 


0.90 


121, 


,51 


± 


0, 


.33 


-0.84 


3 


1/5/99 


51299.62 


5 


26, 


.59 


± 


0, 


.37 


0.59 


121, 


.86 


± 


0, 


.32 


-0.91 


3 


13/5/99 


51311.55 


10 


27, 


.84 


± 


0, 


,14 


0.18 


125, 


.23 


± 


0, 


.20 


0.15 


3 


14/5/99 


51312.55 


10 


28, 


.17 




0, 


,15 


0.43 


124, 


,94 




0, 


.27 


-0.32 


3 


18/5/99 


51316.56 


10 


28, 


.29 




0, 


,14 


0.13 


125, 


,80 




0, 


.15 


-0.18 


4 


28/5/99 


51326.55 


10 


28, 


.83 




0, 


,10 


-0.16 


128, 


.01 




0, 


.13 


0.32 


4 


1/6/99 


51330.52 


10 


29, 


.08 




0, 


,06 


-0.16 


129, 


.27 




0, 


.13 


0.93 


5 


16/6/99 


51345.47 


10 


29, 


.99 


± 


0, 


,12 


0.13 


130, 


,65 




0, 


.16 


-0.08 


4 


21/6/99 


51350.53 


10 


29, 


.48 


± 


0, 


,27 


-0.48 


132, 


,64 




0, 


.28 


1.12 


3 


28/6/99 


51357.48 


10 


30, 


.33 


± 


0, 


.18 


0.31 


133, 


.21 


± 


0, 


.27 


0.62 


3 


28/7/99 


51387.46 


10 


29, 


.36 


± 


0, 


.33 


0.00 


137, 


.39 


± 


0, 


.40 


0.08 


3 


29/7/99 


51388.43 


10 


29, 


.14 


± 


0, 


,14 


-0.17 


137, 


,44 


± 


0, 


.15 


-0.03 


4 


10/8/99 


51400.43 


10 


28, 


.14 


± 


0, 


,30 


-0.52 


139, 


,87 




0, 


.32 


0.41 


3 


28/1/00 


51571.71 


40 


8, 


.81 


± 


0, 


,05 


-0.11 


235, 


.83 




0, 


.12 


-0.89 


3 


2/2/00 


51576.68 


40 


8, 


.60 


± 


0, 


,09 


-0.12 


245, 


.10 




0, 


.12 


-0.55 


3 


6/2/00 


51580.69 


40 


8, 


.26 




0, 


,25 


-0.29 


252, 


,60 




0, 


.25 


-0.51 


2 


16/2/00 


51590.71 


40 


7, 


.20 


± 


0, 


.09 


-0.34 


273, 


,83 




0, 


.20 


-0.37 


3 


18/2/00 


51592.69 


40 


6, 


.82 


± 


0, 


.04 


-0.20 


278, 


.98 


± 


0, 


.16 


-0.40 


5 


2/3/00 


51605.63 


40 


5, 


.67 


± 


0, 


.08 


-0.16 


91, 


.51 


± 


0, 


.31 


-0.16 


5 


14/4/02 


52378.66 


10 


26, 


.57 


± 


0, 


,15 


-0.58 


124, 


,53 




0, 


.15 


0.19 


4 



Based on these estimates, the correlation for either of the 
components is >0.95 for all the baselines employed except 
for 40 m where the correlation is >0.8. For some 67% of the 
observations it follows that Fi and F2 differ by less than 
0.2% and for a further 19% they differ by less than 0.35%. 
Only for 40 m baseline observations do they differ by up to 
1.4%. These estimates justify the assumption that Fi = r2 
for the purpose of determining p and 9. When the system 
is on the meridian the angle ip in equation Q is equal to 
9 since the SUSI baselines are oriented North-South. In the 
analysis the variation of h and 1/) with hour angle was taken 
into account as was orbital motion once a preliminary orbit 
had been established for the system. Initially (3 was included 



as a free parameter but after its value had been established 
with additional observations as discussed in Section [3. 21 the 
fitting procedure was repeated for each night with (3 fixed 
at the value of 0.868. 

The analysis is complicated by the fact that the ob- 
served correlation, and consequently the parameter 7, varies 
during the night due to spatial seeing effects. Depending on 
the seeing conditions the data were corrected for seeing in 
two different ways. 

In SUSI the spatial scale of the wavefront distor- 
tion ro is measured simultaneously with the correlation 
JDavis et al. 199 9a I. On some nights it was possible to con- 
struct calibration curves of correlation versus d/ro, where 



4 J. Davis et al. 



d is the aperture diameter of the interferometer, using val- 
ues near maxima in the correlation to obtain estimates of 
7. These values were corrected for the cosine term in equa- 
tion Q, using prehminary estimates for p and 9, and the 
corrected values were plotted against d/ro. The data were 
fitted using linear least squares (an approximate linear re- 
lationship is predicted by the standard theory of the effects 
of atmospheric turbulence on an interferometer with wave- 
front tip-tilt correction) and the regression line was used to 
calibrate all the data for the night. The corrected data were 
then fitted using equation Q and Fig. shows the results 
for four nights analysed in this manner in which the least 
squares fitted curve is shown with the corrected data. The er- 
ror bars shown in Fig.^are not the individually determined 
values for the uncertainties in the data points since these 
showed large variations but were clearly underestimates of 
the true uncertainty. The error bars correspond to ±10% in 
the correlation value which is believed to be the best esti- 
mate for the average uncertainty in each point. For the figure 
the maximum correlation has been arbitrarily normalised to 
unity. 

Inspection of Fig.^shows that the plots for 10 May 1997 
and 23 April 1999, both obtained with a baseline of 5 m, are 
almost identical. This is to be expected since they are sep- 
arated in time by 1.9973±0.0004 orbital periods (the value 
adopted for the orbital period is discussed in Section [SJ. 

On some nights the correction procedure described 
above could not be used, generally because there was an 
inadequate number of suitable data points to construct a 
calibration curve. In many of these cases the data showed a 
decrease in the measured correlation with time during the 
night, at least in part due to the increase in seeing effects 
with hour angle. In order to fit equation Q to these nights 
a linear decrease in correlation with time was included as a 
free parameter in the fit. In all cases the linear decrease in 
the raw correlation with time was small and <0.014hr~^. 
A more sophisticated model for the instrumental and see- 
ing related trends in the correlation is not justified for the 
present data. While the inclusion of the linear decrease in 
the fit significantly improved the representation of the data, 
it did have a small effect on the values obtained for p and 6. 
The mean difference between the fits, with and without the 
linear term, is 0.29 mas in p and 0.20° in 6. These differences 
are all of the order of, or less than, the formal uncertainties 
in the values from the fits and do not affect the conclusions 
that we draw from the results. 

Tabled lists the values of p and 6 determined for each 
night of observation. There is a 180° ambiguity in the de- 
termination of 6 since the fringe phase can not be deter- 
mined with a two aperture optical interferometer but the 
values listed in TableQare consistent. The uncertainties are 
one standard deviation values determined using the boot- 
strap resampling method (see for example Babu & Feigelson 
1199t)l l. Due to the long period of the orbit (of the order of 
a year) neither p nor 6 change significantly during most of 
the nights on which observations were made. However, an 
iterative procedure was adopted. A preliminary orbital de- 
termination was used to provide the mean rates of change of 
p and 6 for each night and these were then included in the 
fitting procedure. In all cases the rates of change of p and 9 
did not change significantly during the course of the night. 
The mean rates of change were therefore used to obtain a 



single set of values for p and 9, as listed in Table for the 
mean time of observation for each night. 

The orbital fitting program, to be discussed in Section|3 
accepts the position angles and separations and allows for 
relative weights to be assigned to each observational data 
pair. Inspection of the plots of correlation versus hour angle 
for individual nights reveals that the formal errors listed in 
Tableware, in many cases, not a good guide to the weight 
of an individual vector pair. Some nights with few observa- 
tional points have small formal errors whereas other nights, 
with a larger number of observations, have large formal er- 
rors. Weights of 2, 3, 4 or 5 were therefore assigned to each 
night based, not only on the quality of fit of equation Q, 
but also taking into account the quality and number of mea- 
sured values of correlation. The assigned weights are listed 
in Table □ 

3.1 The Visual Companion 

The effect of the presence of the fainter visual companion on 
the observations has been considered. As far as is known this 
third component has shown no significant motion relative 
to the primary pair. It will modulate the observed correla- 
tion, as do the two bright components of the spectroscopic 
binary, but the modulation will be small. This is partly be- 
cause of the magnitude difference but also because of the 
angular separation of the companion from /3 Cen. The an- 
gular separation moves the matched path condition for the 
companion away from that for f) Cen. Thus, for the /3 Cen 
matched path position the correlation for the companion is 
reduced by an amount that increases with baseline. An in- 
vestigation based on the Hipparcos values for the separation 
and position angle of 0'.'9 and 234° for the faint companion 
{Ferryman et al. 1997| l has shown that the modulation will 
only be significant at the shortest baseline of 5 m and even 
there it will not significantly affect the determination of the 
separation and position angle of (3 Cen. 

3.2 The Brightness Ratio 

The ratio of the minimum correlation Cmin to the maximum 
correlation Cmax at a given baseline is given by 

Cmin _ f Fl — (3V2 \ 
Cmax yFl + PV2 J 

Therefore, by measuring the minimum to maximum cor- 
relation ratio, the brightness ratio /3 can be determined if Fi 
and F2 are known. However, if Fi can be assumed equal to 
F2, equation @ simplifies to 

Cmin _ I 1 ~ /3 \ /„^ 
Cma^ ^ ' 

Based on the estimates of the angular diameters of the 
component stars, the corresponding ratio of Fi to F2 at 
the longest baseline (15 m) used in the determination of 
Cmin/Cmax, is >0.998. Equatiou © can therefore be used 
to determine Observations at a baseline of 5 m were not 
included in the determination of (i because of the potential 
effects of modulation due to the fainter companion on the 
determination of Cmin- At longer baselines the modulation 
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Figure 1. Normalised correlation versus hour angle for /3 Cen on (a) May 10, 1997 at a baseline of 5 m, (b) June 9, 1998 at 10 m, (c) 
April 11, 1999 at 15 m, and (d) April 23, 1999 at 5 m. In each case the filled circles represent the measurements and the line is the fit of 
equation to the measurements. 



is greatly reduced, as discussed in the previous section, and 
observations are not significantly affected by it. 

It was not possible to obtain a reliable value for 
Cmin/Cmax usiug the obscrvatious obtained in the manner 
described in Section|5| This was because, for very small val- 
ues of correlation, the on-line correction for the temporal 
effects of seeing fails. A feature of the SUSI data handling 
system is that it provides measurements of the correlation 
for a 1ms sample time (Ci), as well as the correlation cor- 
rected for temporal seeing effects (C). For small values of 
correlation Ci is more reliable than C. Since (i depends crit- 
ically on the correlation ratio, it was decided to base the 
determination of Cmin/Cmax on measurements of Ci and a 
modified observing approach was adopted. On three nights, 
during normal observations of /3 Cen, as minima in corre- 
lation were approached, the OPLC was set to track the es- 
tablished matched path condition through the minima and 
correlation measurements were taken continuously instead 
of measuring delay curves. A quadratic fit to the observed 
values of Ci through each minimum was made to determine 
the minimum value of Ci . To obtain the maximum value of 
Ci, the maximum value of C given by the fit of equation Q 
for each night was scaled by the mean ratio of Ci/C de- 



termined from measurements of correlation away from the 
minima. The ratio Cmin/Cmax was then obtained from the 
ratio of the values of Ci for the mimima and maxima. The 
weighted mean value for the ratio, based on the results for 
5 minima observed over three nights at baselines of 10 m 
and 15 m in April 1999, is 0.0050 ± 0.0012. Substitution in 
equation @ leads to a value for (3 — 0.868 ±0.015. This cor- 
responds to a magnitude difference between the components 
of 0.15 ± 0.02 at the observing wavelength of 442 nm. This 
is in agreement with the original interferometric discovery 
that /3 Cen is a binary with components of approximately 
equal brightness ( |Hanbury Brown, Davis fc Allen igfU l and 
is consistent with the limit on the magnitude difference of 
< 0.3 at 486 nm reported by Robertson et al. 1199911 . 



Fig. |5| shows the values of correlation (Ci) observed 
through a minimum on 12 April, 1999 at a baseline of 15 m. 
The correlation measures have been normalised by the value 
of Ci for the maxima on the night as explained above. 
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12 April 1999: 15 metres 




Hour Angle 

Figure 2. Normalised correlation at a baseline of 15 m versus 
hour angle on April 12, 1999 through a minimum in the correla- 
tion. The filled circles represent the measurements and the line is 
a quadratic fit to the measurements. Further details are given in 
the text. 



4 ANALYSIS OF THE SPECTROSCOPIC 
DATA 

The common parameters determined in independent orbital 
fits to the interferometric and spectroscopic data are tlie or- 
bital period P, epocii of periastron passage T, eccentricity 
e and tiie longitude of periastron uo. An initial orbital fit to 
the interferometric data revealed small but significant differ- 
ences in the eccentricity and longitude of periastron from the 
spectroscopic values published by Ausseloos et al. (120021 . As 
the value of the eccentricity is of major importance in the 
derivation of accurate masses and given that the error esti- 
mates resulting from the VCURVE code used by Ausseloos 
et al. II2UU2II turn out to be too optimistic (Harmanec, Uyt- 
terhoeven & Aerts 2004, Uytterhoeven et al. 2004) an at- 
tempt was made to resolve the small discrepancies by re- 
determining the orbital parameters from the spectroscopic 
data by using a different code. 

Ausseloos et al. 1120021 derived estimates for the orbital 
elements of the primary component relative to the mass cen- 
tre by applying the Lehmann-Filhes method flSSi^ on the 
averaged radial velocities which were determined from two 
Si III triplet lines at 4553A and 4568A. Here we apply the 
publicly available code fotel IIHadrava 19901 on the same 
dataset. The results are listed in the upper part of Table |5| 
and are to be compared with those given in Table 3 of Aus- 
seloos et al. 120021 . As we will see in Section |21 the new 
spectroscopic value for the eccentricity (e = 0.819 ± 0.002) 
is in agreement with the interferometric value. The error es- 
timates by FOTEL, which are based on a covariance matrix 
and should therefore be very reliable, are slightly difi'erent 
from the ones stated by Ausseloos et al. 1201)21 . 

In order to derive the most accurate value for the mass 
ratio, we also redetermined the orbital elements using the 
radial velocities of both the primary and secondary, again 
by applying fotel on the same dataset as used by Ausseloos 
et al. 12UU21 . These results, which are listed in the lower part 
of Table |51 only differ by less than the uncertainties in the 
values of the parameters from the ones obtained previously 



Table 2. The orbital parameters of /3Cen derived from spec- 
troscopy with FOTEL (Hadrava 1990). 



■Porb (days) 


357.00 ± 0.07 


v-y (km s"-"^) 


10.4 ± 3.1 


Ki (km s-i) 


63.8 ± 0.6 


Eo (HJD) 


2451600.39 ± 0.15 


e 


0.819 ± 0.002 


u) (degrees) 


63.6 ± 3.0 


K2 (km s-i) 


63.8 ± 0.8 


Ml /Ma 


1.00 ± 0.04 



with the Lehmann-Filhes method. A remarkable result is the 
equal amplitudes of both orbits, which implies equal masses. 



5 THE ORBITAL PARAMETERS 

The orbital parameters for the /3 Cen system were ini- 
tially determined from the vector separations listed in Ta- 
ble fusing a modification of the program developed at the 
Center for High Angular Resolution Astronomy (CHARA) 
jHartkopf, McAlister fc f 'Vanz 19891 . The modified program 
uses the simplex algorithm ( Ncldcr & Mead 1965 ) instead of 
a grid search to determine the best fitting orbit and deter- 
mines the uncertainties in the orbital elements using Monte 
Carlo simulation. Both programs accept the position angles 
and separations and allow for relative weights to be assigned 
to each observational data pair. The programs include cor- 
rections for precession and proper motion. As discussed in 
Section|3 weights of 2, 3, 4 or 5 were assigned to each night 
of observation based on the quality of the fit of equation Q 
but also taking into account the quality and number of mea- 
sured values of correlation. The assigned weights are listed 
in Table Both programs were run with the mean MJDs in 
Table converted to Besselian epochs. There are no signifi- 
cant differences between the orbital elements determined by 
the two methods. 

The orbital parameters obtained from the fit to the 43 
entries in Table using the modified program, are listed 
in Table |3 The common parameters determined spectro- 
scopically and listed in Table |21 are repeated in Table |3| 
The common parameters are in good agreement but we note 
that the orbital period from the spectroscopic solution has 
a much smaller uncertainty than that determined interfer- 
ometrically. This is to be expected since the spectroscopic 
observations spanned 12 years compared with ~ 5 years 
for the interferometric observations. Following the example 
set by Herbison-Evans et al. 119711 in their analysis of in- 
terferometric observations of a Vir, the orbital period was 
fixed at the spectroscopic value and the orbit fitting pro- 
gram was run again on the interferometric data. The re- 
sulting values for the remaining six parameters are listed in 
Table 13 The differences between the parameters, with and 
without the period fixed, are generally less than, or of the 
order of the uncertainties in the values of the parameters. 
The fit with the period fixed is marginally better than with 
the period a free parameter and, although it is not evident 
in the rounded values shown in Table El the uncertainties 
in the derived parameters were all smaller with the period 
fixed. For this reason and the longer span of the spectro- 
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Table 3. Orbital parameters for (3 Cen. The first interferometric column is for a fit for all seven parameters; 
the second interferometric column is for a fit with the orbital period fixed at the spectroscopic value; the 
spectroscopic values are from lAusseloos et al. 20021 . Further details are given in the text. 



Parameter 


Value 


Value 


Value 




(Interferometric) 


(Spectroscopic) 


(Interferometric) 








(Spectroscopic P) 


Period (P) 


35r?40±0'?20 


357?00±0'?07 


357? 00 


Angular semi-major axis (a) 


0"02509±0^'00025 parameters 




0y02532±0^'00023 


Inclination (i) 


67?1±0?5 




67?4±0?3 


Epoch of periastron passage (T) 


2451599'? 9±0'?3(HJD) 


2451600? 39±0?15 (HJD) 


2451600? 0±0?3 (HJD) 


Position angle of line of nodes (Q) 


288?5±0?4 




288?3±0?4 


Eccentricity (e) 


0.821±0.005 


0.819±0.002 


0.824±0.004 


Longitude of periastron (uj) 


60?8±0?6 


63.6±3.0 


61?3±0?6 



scopic observations, the orbital parameters for the fit with 
the orbital period fixed are adopted as best representing the 
interferometric data. 

The observationally determined orbital points and the 
orbit for the adopted orbital parameters are plotted in Fig:|3 
following the common double-star convention (IHeintz 197811 . 
Also plotted are the positions on the fitted orbit for the 
epochs of the individual observational points. Error ellipses 
for the observational points have not been plotted since they 
are comparable in size or smaller than the symbols represent- 
ing the observational points and, as discussed in Section |21 
they are not always the best representation of the weight 
of an individual point. We note that the orbit with the or- 
bital period as a free parameter is almost indistinguishable 
from the plotted orbit, except that many of the observational 
points lie a little further from the corresponding positions on 
the orbit. This is further evidence supporting the decision 
to accept the orbital fit with the period fixed. 

In Fig.|21the vector position determined with MAPPIT 
by Robertson et al. II1999II has been plotted and, while it was 
not included in the fit, it is seen to be consistent with the 
orbit derived from the SUSI data. There is a 180° ambiguity 
in the position angle of the MAPPIT point relative to the 
SUSI data and in the figure the point has been plotted at 
position angle 172°rather than the value of 352°given by 
Robertson et al. 

Examination of Fig. |3 reveals that, while over most of 
the orbit the observational points are distributed on either 
side of the curve, the observations made in 2000 near to 
periastron all lie inside the fitted orbit. This is shown on 
an expanded scale in Fig. 01 where a section of the orbit 
close to periastron is plotted on an enlarged scale. While 
we have not been able to provide an explanation for this 
phenomenon we conjecture that, since the separation of the 
component stars is less than 7 stellar radii at periastron, 
the possibility of interaction may affect the brightness dis- 
tributions such that the components appear closer than their 
centres of mass. Miroshnichenko et al. (I2001|l have reported 
interaction around periastron for the early B type binary 5 
Sco involving a Be outburst. Like /3 Cen, 5 Sco has a highly 
eccentric orbit, but with an orbital period of ~ 10.6 years. 
Although 5 Sco spends longer close to periastron, the clos- 
est approach is of the order of 30 times the primary radius 
compared with less than 7 radii for /3 Cen. This suggests 
that some form of interaction is likely to occur. 



-25 




5 ' 
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y (mas) 

Figure 3. The orbit of (3 Centauri. The data points are the SUSI 
results from Table 111 except for the MAPPIT point which is from 
Robertson et al. (1998). SUSI observations: 1997 0; 1998 A; 1999 
O; 2000 •; 2002 □. MAPPIT observation: 1995 X. The points on 
the fitted orbit that correspond to the observational points are 
marked with +, Further details are given in the text. 



6 COMBINATION OF THE SPECTROSCOPIC 
AND INTERFEROMETRIC RESULTS 

The interferometric and spectroscopic results can be com- 
bined to yield the distance to the system, the absolute mag- 
nitudes of the components and the masses of the compo- 
nents. Not all the parameters are required for this purpose 
and the relevant parameters from the two techniques are 
summarised in Table 2] A mean of the interferometric and 
spectroscopic values for the eccentricity is adopted. 

As the two sets of data result in a single self-consistent 
set of parameters we are able to derive the masses of the 
components with high precision. Functions containing the 
masses of the primary (Mi) and secondary (M2) components 
in solar mass units and the semi-major axis (a) of the orbit 
in AU are given by IIHeintz 1978t : 

Ml sin^ i = 1.036 X {Ki + Ki fK-zPil - e'^ f^^, (7) 
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Figure 4. A section of the orbit of /3 Centauri around periastron. 
The data points are the SUSI results from Table [Tl Each point 
has its error ellipse plotted but, in all but one case, it is too small 
to extend outside the plotted point. The curve is the orbital fit 
shown in Fig. |3] The location of periastron is marked by a circle. 

Table 4. Values of parameters for f} Cen determined directly 
by intcrferometry and spectroscopy and used to derive the basic 
physical parameters of the system. 



Parameter Value Source 



Period (P) SST^OOitO^O? S 

Ki (km s-i) 63.8±0.6 S 

K2 (km s-i) 63.8±0.8 S 

Eccentricity (e) 0.821±0.003 I + S 

Angular semi-major axis (a) 0"02532±0^'00023 I 

IncUnation (i) 67?4±0?3 I 

Brightness ratio (/3) (at A442 nm) 0.868±0.015 I 



Ma sin^ i = 1.036 x IQ^\Ki + K2fKxP{l ~ e^f^, (8) 
asin j = 9.1919 x \Q-^{Kx -f K2)P{\ - ef^. (9) 



The values for these functions, computed from the data 
in Table H are listed in Table 1^1 

6.1 The distance to fi Cen 

The parallax of the system can also be determined by com- 
bining the interferometric and spectroscopic results and is 
given by: 



(asini) 

Table 5. Values of parameters for Cen derived from the pa- 
rameters listed in Table HI using equations (l)-(3). M\ and Afa 
are the masses of the primary and secondary respectively and a 
is the semi-major axis of the orbit. 



Parameter 


Value 


Reference 


Misin^i (M©) 


7.15±0.24 


Equ. Q 


Masin^i (Mq) 


7.15±0.23 


Equ. 


asini (AU) 


2.391±0.032 


Equ. © 



A parallax determined in this way has been appropri- 
ately termed an orbital parallax ( [Armstrong et al. 1992||. 
Combining a" and i from Table |1] with asini from Table |3 
gives the orbital parallax of /3 Cen equal to 9.78±0.16 mas. 
This is to be compared with the HIPPARCOS parallax of 
6.21±0.56mas l |Perryman et al. 1997| l. There is a large and 
significant diflerence in spite of the relatively large fractional 
uncertainty in the HIPPARCOS value. The HIPPARCOS 
catalogue recognises [i Cen as a binary system but has no 
mention of the dual nature of the primary component of 
/3 Cen and lists only the faint third magnitude companion. 
We conjecture that the HIPPARCOS value for the parallax 
is seriously affected by the binary nature of the primary of 
/3 Cen. In fact, the HIPPARCOS observations of spectro- 
scopic binaries were not fitted in general with an orbital 
model to derive the value for the parallax. While repro- 
cessing of the HIPPARCOS astrometric data has been per- 
formed for many single-lined spectroscopic binaries (see for 
example (IPourbaix fc Boffin 20031 and references therein), 
it has not been done for /3 Cen. For this reason we adopt the 
value derived directly from the interferometric and spectro- 
scopic data. The parallax we have adopted is equivalent to 
a distance of 102.3±1.7pc. 

6.2 The absolute magnitudes and spectral types 
of the components of /3 Cen 

Breger llTMTl has given the following photometric data for 
/3 Cen: Y = 0.61, (B-V) = -0.24 and (U-B) = -0.99. Since 
the magnitude difference is only 0.15±0.02 at a wavelength 
of 442 nm (from the brightness ratio of 0.868±0.015 given in 
Table [IJ and, as will be shown in Section 16.31 the masses 
of the components are equal, any difference in the colours 
of the two components will be very small and they will be 
assumed to be the same. Reddening will be small for this 
relatively nearby system (distance ~ 102 pc) and E(B-V) is 
estimated to be ~ 0.03. It follows that {B-V)q = -0.27 and 
we adopt Av = 0.09 ± 0.03. 

Using the parallax from Section fO. li the absolute visual 
magnitude Mv of the system is given by: 

Mv = V - 0.09 - 5.0 log(102.4/10). (11) 

Thus Mv = -4.53±0.05 for the two components com- 
bined. The brightness ratio of 0.868±0.015, although mea- 
sured at A442 nm, is unlikely to differ significantly in the 
V band because of the equality of their masses and the 
small difference in luminosity between the two components. 
Assuming that the same brightness ratio is valid in the vi- 
sual band, the absolute visual magnitudes for the two com- 
ponents are Myprlmary ~ —3.85 ± 0.05 and Mvsocondary ~ 

-3.70 ±0.05. 

These values for the absolute visual magnitudes lie be- 
tween the values given by Balona & Crampton 119741 for 
spectral types BlIII (-4.1) and B2III (-3.7) and are close to 
the value of -3.9 given by Schmidt-Kaler (11^5^ for B2III. 
They lie in the middle of the range of absolute visual mag- 
nitudes for /3Cep stars ( |Sterken fc Jerzykiewicz T993J. 

6.3 The masses of the components of (5 Cen 

The combination of the value for i from Table 2] with the 
values for Afi sin^ i and M2 siv? i from Table |3 gives the 
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Table 6. The parallax (vr), masses (Mi & M2) and absolute visual 
magnitudes {Miy & M2v) for f) Cen and its component stars. 
The absolute visual magnitudes were computed directly from the 
photometric data, the adopted parallax, and the brightness ratio 
P given in Table 1^ 



Parameter 


Value 


vr (mas) 


9.77±0.15 


Ml (M©) 


9.09±0.31 


M2 (Mq) 


9.09±0.30 


Miy 


-3.85±0.05 


M2V 


-3.70±0.05 



mass of the primary Mi equal to 9.O9±O.32M0. Similarly 
the mass of the secondary M2 is 9.09±0.31 M©. These are 
the most accurate mass estimates of any /3Cep star in a 
binary system known to date. 

Masses for early type giants are in general rather uncer- 
tain. The mass for a BO III type star is usually taken as ~ 
20 Mq and, for a B5 III type, ~ 7 Mq JSchmidt-Kaler 19821 . 
The values we have determined for the components of /3 Cen 
are significantly lower at BlIII than implied by those gen- 
erally accepted for B giant stars. 

Refined mass estimates at ~3.5 per cent accuracy, 
which is the level of our results for 13 Cen, have recently 
become available for two single (3 Cep stars from astero- 
seismology, i.e. from detailed interpretation of the mea- 
sured frequencies and mode identification of their oscilla- 
tions. Such studies have led to very similar mass estimates 
for these two single stars, i.e. a mass range [9.0,9.5] Mq 
for HP 129929 HAerts et al. 200311 and [9.0,9.9 ] Mq for v Eri 
jPamyatnykh, Handler fc Dziembowski "20041 . The mass of 
the echpsing binary (3 Cep star 16 Lac is also rather well 
constrained and lies in [9.0,9.7] Mq (Tho ul et al. 2 0Q3)l. 

The masses of all the known /3Cep stars as a group of 
variables have mainly been derived from photometric cali- 
brations, e.g. of the Walraven, Ceneva and Stromgren sys- 
tems. This method provided systematically too high masses 
for the three stars with accurate seismic mass determination: 
M ~ 10 to 11.3 Mq for HD 129929; M ~ 10 to 10.8 Mq 
for u'Eri and M = 10 Mq for 16 Lac (values taken from 
Heynderickx, Waelkens & Smeyers 1994; see also Sterken 
& Jerzykiewicz 1993) to be compared with the seismic val- 
ues listed above. The situation is of course worse when we 
consider the mass estimates of photometric calibrations for 
spectroscopic binaries since the binarity is not usually taken 
into account in such global calibrations. Heynderickx et al. 
I|1994ll . for example, list a mass between 14.4 and 15.2 Mq 
from treating (3 Cen as a single star, as they did all other 
stars in their sample. Their estimate for the mass is far too 
high compared with our results. Since a large fraction of the 
/3Cep stars turn out to be a member of a spectroscopic bi- 
nary with not too different components (Aerts & De Cat 
2003), it is unavoidable that the mass estimates of these 
group members from photometry are very unreliable. 



7 SUMMARY 

The orbital parameters for the primary component of (3 
Cen have been determined from interferometric observations 



made with SUSI. These parameters have been combined 
with the high resolution double-lined spectroscopic data for 
/3 Cen reported in Ausseloos et al. II2UU2II . the analysis of 
which has been revised here. The combination of data has 
yielded the masses and absolute visual magnitudes of the 
component stars and an accurate distance to the system. 
The very precise results of this powerful combination of tech- 
niques are summarised in Table |H| 

The results presented in this paper constitute a very 
suitable starting point for any future asteroseismic mod- 
elling of /3Cen's two pulsating components. The star is also 
an ideal target to investigate possible tidal effects on the 
oscillations in view of the very close passage (less than 10 
times the radii of the component stars) of the stars at peri- 
astron although, for this, long-term monitoring over several 
orbits is required. 
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